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■ ABSTRACT 

We investigate the influence of dark energy on structure formation, within five different cos- 
q . mological models, namely a concordance ACDM model, two models with dynamical dark 

5— i ■ energy, viewed as a quintessence scalar field (using a RP and a SUGRA potential form) 

| and two extended quintessence models (EQp and EQn) where the quintessence scalar field 

interacts non-minimally with gravity (scalar-tensor theories). We adopted for all models the 
normalization of the matter power spectrum eg to match the CMB data. In the models with dy- 
namical dark energy and quintessence, we describe the equation of state with wq ps —0.9, still 
within the range allowed by observations. For each model, we have performed hydrodynami- 
\Q | cal simulations in a cosmological box of (300 Mpc ft, -1 ) 3 including baryons and allowing for 

. cooling and star formation. The contemporary presence of evolving dark energy and baryon 

C*~) 1 physics allows us to investigate the interplay between the different background cosmology 

| and the evolution of the luminous matter. Since cluster baryon fraction can be used to con- 

qq . strain other cosmological parameters such as Q m , we also analyse how dark energy influences 

" the baryon content of galaxy clusters. We find that, in models with dynamical dark energy, the 

f^*) , evolving cosmological background leads to different star formation rates and different for- 

mation histories of galaxy clusters, but the baryon physics is not affected in a relevant way. 
We investigate several proxies for the cluster mass function based on X-ray observables like 
temperature, luminosity, M gas , and Yx- We conclude that the X-ray temperature and M gas 
functions are better diagnostic to disentangle the growth of structures among different dark 
energy models. We also evaluate the cosmological volumes needed to distinguish the dark en- 
ergy models here investigated using the cluster number counts (in terms of the mass function 
and the X-ray luminosity and temperature functions). Relaxed, massive clusters, when studied 
in regions sufficiently far from from the centre, are built up in a very similar way despite the 
different dark energy models here considered. We confirm that the overall baryon fraction is 
almost independent of the dark energy models at a few percent level. The same is true for the 
gas fraction. This evidence reinforces the use of galaxy clusters as cosmological probe of the 
matter and energy content of the Universe. 
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1 INTRODUCTION 



Over the last decade great observational evide nce ( Riess et al 



19981 : IPerlmutter etafl 1 19991 : Ijarosik et alj 1201 it IVikhlinin et al 
2009) has shown that at the present time the Universe is expand 
ing at an accelerated rate. This fact can be attributed to a compo 
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nent with negative pressure, which is usually referred to as dark 
energy, that today accounts for about 3/4 of the entire energy bud- 
get of the Universe. The simplest form of dark energy is a cos- 
mological constant term A in Einstein's equation, within the so- 
called ACDM cosmologies. Though in good agreement with ob- 
servations, a cosmological constant is theoretically difficult to un- 
derstand in view of the fine-tuning and coincidence problems. A 
valid alternative consists in a dynamical dark energy contribution 



2 



that changes in time and space, often associated to a scalar field 
(the 'cosmon' o r 'quintessence') evolv ing in a suitable potential 
dWetteric r] |l988l : lRatra& Peebles! 1 1988h . Dynamical dark energy 
allows for appealing scenarios in which the scalar field is the medi- 
ator of a fifth for ce, either within scalar-tensor theories or in inter- 
acting scenarios dWetterichll 1 9951; I Amendolall20qd: I Boisseau et al.l 



2000; Pettorino & Ba ccigalupil 12008; Mota et al. 2008], and refer- 



ences therein). In view of future observations, it is of fundamen- 
tal interest to investigate whether dark energy leaves some im- 
prints in structure formation, giving a practical way to distin- 
guish among different cosmologies, as recently investig a ted in 
iHu & Sawickj j2007h.lB~aldi et all J20I0I) . IZhao et aTl j2010hjBaldil 
j201 lh . lBaldi & Pettorinol feOl lh and lwintergerst et alj d2010h . ~ 

In this paper, we study the general properties of galaxy clus- 
ters in different dark energy cosmologies. Galaxy clusters are the 
largest virialized objects in the Universe and are considered to be 
a fair sample of the overall matter distribution of the Universe it- 
self. They contain a large amount of gas in the form of diffused 
ionized plasma known as intracluster medium (ICM), which emits 
in the X-ray band. The X-ray properties of galaxy clusters such as 
luminosity and temperature trace the total mass of the cluster itself, 
and hence can be used to study global properties of these objects. A 
lot of observational work (Chandra, XMM-Newton) has been made 
in recent years, and future missions (e.g. IXO, eROSITA, WFXT) 
are under study to improve the characterization of these objects 
in the X-rays. The properties of galaxy clusters, in particular their 
mass, can be investigated also in the optical region of the spectrum 
through gravitational lensing, which gives independent estimates 
from X-rays. 

In this paper, we analyse the properties of simulated galaxy 
clusters in cosmologies with different dark energy models. We fol- 
low the formation and evolution of structures in hydrodynamical 
simulations of a cosmological box of size (300 Mpc /i -1 ) 3 for five 
different cosmologies. These cosmologies are in general character- 
ized by the presence of a dynamical dark energy component, i.e. a 
dark energy component with density and equation of state evolv- 
ing with time, and they will be introduced and discussed in detail 
in Sect. [2] After presenting the selection and composition of the 
sample we extract from the simulations, we will study the mass 
function at different redshifts. Different dark energy models have a 
different CMB normalization of the spectrum of the perturbations 
and a different growth factor, both things affecting the mass func- 
tion of galaxy clusters. Then we will study the X-ray luminosity 
and temperature functions, because both quantities are a proxy for 
the cluster mass and can therefore be used to search for dark en- 
ergy imprints. In this case the advantage is that these quantities are 
somehow directly observable with space facilities, while evaluating 
the mass of a galaxy cluster requires the assumption of a model or 
physical hypotheses (such as hydrodynamical equilibrium) that can 
introduce systema tics (for example i f the cluster is not in a relaxed 
state), as shown in iRasia et al.l (2006). We will also examine the X- 
ray luminosity-temperature (L — T) relation of our sample, in order 
to check whether there is any clear discrepancy between the proper- 
ties of the simulated objects in comparison with the observe d ones 
dArnaud & Evrardlll999i : lMaughanll2007l ; IPratt et alj|2009l) . After 
the analysis of the global properties of the simulated sample, we 
will concentrate on the internal properties of the single clusters and 
study the relative distribution of the different mass components. In 
particular we will check the dark energy dependence of the baryon 
fraction fb ar - The study of the baryon fraction in galaxy clusters, 
either in the form of ICM or of stars in galaxies, is useful to un- 
derstand the formation history and the properties of these objects. 



Combined with other independent information, such as for exam- 
ple the value of the baryonic density Qb derived from primordial 
nucleosynthesis, ft, ar can be used to derive the value of the matter 
density of the Universe, fl m — Qbfbar- 

The paper is organized as follows. In Sect.[2]the different cos- 
mological models will be introduced; in Sect. [3] the simulation set 
will be discussed and in Sect. [4] the study of the mass functions 
of the selected sample will be analysed. Sect. [5] is centred on the 
analysis of the X-ray L — T relation, while in Sect.|6]the X-ray ob- 
servables functions will be studied. In Sect. [7] we will describe the 
analysis and the results of the study of the cluster baryon fraction, 
while conclusions will be drawn in Sect. [8] 



2 THE COSMOLOGICAL MODELS 

In a homogeneous and isotropic universe, the cosmological expan- 
sion depends on the various energy components, as described by 
the Friedmann equation: 



H 1 



Hi 



floi exp 



1 + Wi(a') 1 
H — -da 



+ 



(i- 



(1) 



Here, the Hubble parameter is defined as H = a /a (where the 
dot denotes the derivative with respect to the cosmic time t), 
Qoi = poi/poc is the current density parameter of the i-th com- 
ponent of the universe, po c = 3Hq / 8nG is the critical density at 
the present time, Wi is the equation of state parameter of the i-th 
component (u>i = Pi/pi) and the sum is taken over all compo- 
nents. We have used natural units and set c = 1. The first term in 
equation |QJ includes densities associated to each constituent of the 
universe while the second term accounts for any possible deviation 
from flat geometry. We can assume that the universe is constituted 
of three different components: matter [baryons and cold dark mat- 
ter (CDM)], for which w m = 0; radiation (photons plus relativistic 
matter) with w r = 1/3, whose contribution is nowadays negligi- 
ble; dynamical dark energy (DE), which, in the simplest case, be- 
haves as a fluid with negative time dependent ibde and provides 
the present accelerated expansion of the universe. We consider flat 
cosmological models for which ^ - Ooj = 1, so the curvature term 
in the Friedmann equation will be equal to zero. With this in mind, 
and introducing n = 8nG, we can express equation Q} as 



H 1 



K \ ^ 

3 ^ 



poi exp —3 



1 + Wi(a') 



da 



(2) 



For our analysis we have considered three possible sets of cos- 
mological models. The first is the standard ACDM model, that we 
use as a reference model, where dark energy is represented by the 
cosmological constant. This model is in agreement with present ob- 
servations, though theoretically it is intrinsically affected by fine- 
tuning and coincidence problems. Alternatively, dark energy could 
be a dy namical component, seen as a scala r field rolling down a po- 
tential jWetterich|[l988l ; lRatra & Peeblesll 19881) . If the scalar field 
is minimally coupled to gravity, this class of scenarios is still af- 
fected by fine-tuning and coincidence problems, as much as in the 
standard ACDM model. However it is interesting, for our analysis, 
to consider such dynamical cases, where a time varying equation 
of state is present. Numerical simulations of quin tessential cold 
dark matter have been presented, for example, in Ijennings et al.l 
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bold) . More interestingly, the dynami cal scalar field c o uld be cou- 
pled to oth er species, as addressed inlWetterichl J 1 995t) , lAmendolal 
(2000) and lPettorino & Baccigalupil d2008h . We limit ourselves to 
the case in which the coupling inv olves universally all s pecies, as it 
happens in scalar-tensor theories dBoisseau et alj|20po|). The latter 
have b een also investigated with in F(R) theories in lSchmidt et al.l 
d2009h and lOvaizuetai] (§008). Af-body simulations with a cou- 
pling to g ravity (extend ed quintessence) have been studied, for ex- 
ample, in iLi et alj l201lh . Note that hydro- sim ulations including a 
coupling to dark matter have been presented in lBaldi et alj d2010l) . 
iV-body simulations for cou pled dark energy hav e b een investi- 
ated in lMaccio et al.1 l2004h and lZhao et ailfcoiol) . In lLi & Zhad 
> et al.1 



2009) and Zhao et al. (2010) the effect of scalar field perturbations 



was also take n into account. Fu rther fifth force couplings have been 
simulated in iLi & Zhaol d2009l) . The impact of early dark energy 
in stru cture formation has been considered in iGrossi & Springell 
(2009). Our analysis differs from previous works because of dif- 
ferent dynamics in the dark energy models considered and of the 
contemporary inclusion of baryonic physics in the framework of 
evolving dark energy scenarios. 

We will now specify in more detail the cosmologies here con- 
sidered and our choice of the parameters. 



where V(<p) is the potential in which the scalar field (f> rolls. 

Similarly, the evolution of the quintessence scalar field <j) is 
given by the Klein-Gordon equation, obtained from equation (O 
substituting equations Q and {8j, 



; + 3^ + ^ = o. 



(9) 



We note from equations 10 and ([8} that when the kinetic term </> 2 /2 
is negligible with respect to the potential term V ((/>), then w^, — > 
— 1 and the ACDM case is recovered. 

In this paper, as potentials for minimally coupled quintessence 
models, we consider an inverse power law potential 



v(4>) 



M 4+a 



(10) 



the so called RP potential dRatra & Peebles! Il988h , as well 
as its generalizatio n suggested by supergravity arguments 
jBrax & MartirJI 19991) . known as SUGRA potential, given by 



M 4+a 

V W = exp(4^G0 2 ) , 

where in both cases M and a > are free parameters. 



(11) 



2.1 ACDM 

The first model considered is the concordance ACDM model, with 
the values of the cosmological parameters taken from WMAP3 (see 
Subsect l2.4l (. This model is characterized by the presence of a dark 
energy component given by a cosmological constant A, with a con- 
stant wa — —1, so that equation (|2} can be written as 



H 1 



P0r, 



(?)' 



+ POr 



+ P0A 



(3) 



with m, r, and A denoting the matter, radiation and cosmological 
constant components respectively. 



2.2 Minimally coupled quintessence 

The second case that we consider here is that of a dynamical dark 
energy, given by a quintessence scalar field <t> with an eq uation of 
state w = w(a) dWettericrJ[l988l ; iRatra & Peebleslll988n . In gen- 
eral, we can express equation ^ as 



/'(in 



a \ 3 , /a<^ 4 
— + POr- — 
a 1 \ a 



+ P<j> 



(4) 



where (f> denotes the quintessence component. The evolution of p^ 
can be obtained from the continuity equation 



P4, + iH(p 4> + p$) = , 
so that formally 

,, ra 1 + W4,(a') 

P4 = Po0 exp 



da 



(5) 



(6) 



where = / p$. The conserved energy density and pressure of 
the quintessence scalar field that appear in equation lO are defined 

as 



= ^ +v(<t>) 



(7) 
(8) 



2.3 Scalar-tensor theories 

It is well possible that the quintessence scalar field might interact 
with other species. Here we consider the case in which <j> interacts 
non-minimally with gravity dWettericbll988l ; lBoisseau et al.ll2000h 
and we refer in particular to th e exten d ed quintessence (EQ ) mod- 
els described in | Perrotta et alj feOQO), Pet torino et alj d2005h and 
|Pettorino~& Baccigalupil l2008h . Scalar-tensor theories of gravity 
are generally described by the action 



2k 



+ £ fluid , (12) 



where R is the Ricci scalar, the function f(<f), R) specifies the cou- 
pling between the quintessence scalar field and the Ricci scalar, 
(jj(<f>) and V(4>) specify the kinetic and potential terms respec- 
tively and the Lagrangian Cauid includes all the components but 
(f>. Here we assume for the sake of simplicity a standard form for 
the kinetic part, uj(4>) = 1, and we define the coupling function as 
f(4>, R) — K,F(<j})R, where F(<j)) is chosen to be 



(13) 



Here k = 8ttG*, where G* represents th e "bare" gravitational con- 
stant dEsposito-Farese &Polarskill200lh . which is in general dif- 
ferent from the Newtonian constant G and is set in such a way 
that locally 1/k + £(<^> 2 — 2 ) = 1/8ttG in order to match local 
constraints on General Relativity. The parameter £ represents the 
"strength" of the coupling. In particular we consider here a model 
with positive coupling £ > (EQp) and one with negative £ < 
(EQn). Note that a dependence on the sign is expected in this case 
dPettorino~& Bacciaalupi 2008). The limit of General Relativity is 
recovered when cj.jbd 2> 1, where 



W.7BD 



F(4) 



[dFW/d<t>\ 



(14) 
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String ent constraints for this quantity come from the Cassini mis- 
sion l lBertotti et alj|2003h on Solar System scales, where wjbdo > 
4 x 10 4 . However, it has been noted that such con straints may not 
apply at cosmological scales dClifton et al J 120051) where comple- 
mentary bounds, obtained combining WMAP1 and 2dF large scale 
structure data, pr ovide the less tight lim it of cjjbdo > 120 at 95% 
confidence level dAcquaviva et alj|2005l) . 

In EQ models, we can d efine a conserved density and pressure 
for the scalar field, given by dPerrotta & Bacciga lupi 2002): 



1 ■ 



<f + V(<j>) - 3HFty) + 3iT 



F{<)>) 



P4> = 2^ ~ V(4>) + F(<j>) + 2HF{4>) + 



(2H + 3H 2 ) 



■F(4>) 



(15) 



(16) 



respectively. The evolution of the scalar field follows from the 
Klein-Gordon equation 



■ + Wtj) + 



dV{<j>) 1 dF(<j>) 



R. 



(17) 



where the Ricci scalar is given by 

R = 6(H + 2H 2 ). (18) 
As in previous scenarios, equation (|2j can be again expressed as 



H 2 = K - 



par, 



(S)' +*(?)' 



+ p<l> 



(19) 



where is the conserved energy density defined in equation H5\ . 
In this paper, as underlying potential for the extended quintessence 
models, we use the RP potential in equation l |10| >. Looking at equa- 
tions j 1 3b . d!5l >, j 1 6b we notice that minimally coupled quintessence 
is recovered for £ — ► 0. 

For an extensive linear tre atment of EQ models we refer to 
IPettorino & Baccigalupil J2008). Here we only recall for conve- 
nience that EQ models behave like minimally coupled quintessence 
theories in which, however, a time dependent effective gravitational 
interaction is present. In particular, in the Newtonian limit, the Eu- 
ler equation for CDM can be written as 



Vt>m + HVv 



iwGM m S(0) 



. 



(20) 



in terms of the cosmic time, where we have redefined the gravita 
tional parameter as 

2[F + 2(dF/dcf>) 2 



G = 



1 1 



(21) 



[2F + 3(dF/d<j>) 2 } 8ttF ' 

The latter formalism is general for any choice of F(<j>). For the 
coupling here chosen and given by equation J 1 3 1 > we have 

+ (l + 8£)£0 2 - C^o] 1 _ 



G = 



i 

8irG« 



8tvG, 



+ (1 + ~ €08 + 8<(0 2 



For small values of the coupling, that is to say ( < 1, the latter 
expression becomes 



(23) 



which manifestly depends on the sign of the coupling £. We note 
that, since the derivative of the RP potential in equation l |10t with 
respect to cj> is dV(4>)/d<f> < 0, we have (f> 2 < <j>^. This leads to the 
behaviour of G/G, shown in Fig.Q] 




Figure 1. Correction to the gravity constant for the two extended 
quintessence models, EQp (cyan) and EQn (red), as expressed in equation 
423t . Note that the corrections are only within the percent level. 



2.4 Choice of the parameters 

As a reference model we use the A CDM model, adapted to the 
WMAP3 values dSpergel et ai]|2007ri . with the following cosmo- 
logical parameters: 

• matter density: Qom = 0.268 

• dark energy density: f2oA = 0.732 

• baryon density: fiot = 0.044 

• Hubble parameter: h — 0.704 

• power spectrum normalization: ag = 0.776 

• spectral index: n s = 0.947 

We trimmed the parameters of the four dynamical dark energy 
models so that too = w(Q) w —0.9 is the highest value still con- 
sistent with observational constraints in order to amplify the effects 
of dark energy. Fig. [2] shows the evolution with redshift of w in 
each cosmology. The parameters f2 Um , f^oA, fiof» h, and n s are the 
same for all the models, but since we normalize all the dark energy 
models to CMB data from WMAP3, this leads to different values 
of as for the different cosmologies: 



D 4 



C8,DE = 0"8,ACDM - 



i(zcmb) 



D 



+ ,DE 



(«cmb) 



(24) 



assuming zcmb = 1089. This fact, along with the different evo- 
lution of the growth factor D+ (shown in Fig. [3), has an impact 
on structure formation. Table Q] lists the parameters chosen for the 
different cosmological models. 



3 NUMERICAL SIMULATIONS 

In order to study the formation and evolution of large scale 
structures in these different cosmological scenarios we use iV- 
body + hydrodynamical s imulations don e with the GADGET-3 
code dSpringel et al.l l200ll : ISpringell [2005), which makes use of 
the en tropy-conserving formulation of SPH dSpringel"& Hernquist 
2002). We extended t he dark energy implementation as described 
in iDolag et alj d2004l) to allow the code to use an external, tabu- 
lated Hubble function as well as a tabulated correction to the gravity 
constant needed for the extended quintessence models. The hydro- 
dynamical simulations include radiative cooling, heating by a uni- 
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Table 1. Parameters chosen for the different cosmological models: a is the 
exponent of the inverse power law potential; £ is the coupling in the ex- 
tended quintessence models; wjbdo i s the present value of the parameter 
introduced in equation 114) ; wq is the present value of the equation of state 
parameter for dark energy; erg is the normalization of the power spectrum 
as in equation 124) . 



Model 


a 




■fJBDO 


w 


CT8 


ACDM 








-1.0 


0.776 


RP 


0.347 






-0.9 


0.746 


SUGRA 


2.259 






-0.9 


0.686 


EQp 


0.229 


+0.085 


120 


-0.9 


0.748 


EQn 


0.435 


-0.072 


120 


-0.9 


0.729 




Figure 2. Redshift evolution of the equation of state parameter w for 
the different cosmological models considered: ACDM (black), RP (blue), 
SUGRA (green), EQp (cyan), and EQn (red). 



o 
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Figure 3. Redshift evolution of the growth factor D+, normalized to 
the corresponding trg, for the different cosmological models considered: 
ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn (red). 



form redshift-dependent UV background faaardt&Madaull 19961) . 
and a treatment of star formation and feedback processes. The pre- 
scription of star formation we use is based on a sub-resolution 
model to account for the multi-phase structure of the interstellar 
medium (ISM), where the cold phase of the ISM is the reservoir 
of star formation ( Sprin gel & Hernquistll2003l) . Supernovae (SNe) 
heat the hot phase of the ISM and provide energy for evaporat- 
ing some of the cold clouds, thereby leading to self-regulation of 
the star formation and an effective equation of state to describe 
its dyna mics. As a phenomen ological extension of this feedback 
scheme, ISpringel & Hemq uist ( 2003) also included a simple model 
for galactic winds, whose velocity, v w , scales with the fraction r\ 
of the Type II SN feedback energy that contributes to the winds. 
The total energy provided by Type II SN is computed by assuming 
that the y are due to e xploding massive stars with mass > 8M© 
from a Salpete \ dl955h initial mass function, with each SN re- 
leasing 10 51 ergs of energy. We have assumed r\ — 0.5, yielding 
v w ~ 340kms _1 . 

We simulated a cosmological box of size (300 Mpc /i -1 ) 3 , 
resolved with (768) 3 dark matter particles with a mass of itldm ~ 
3.7 x 1O 9 M0 h' 1 and the same amount of gas particles, havin g 
amass of m gas « 7.3 x 1O 8 M hT 1 . As in bolag et ail d2004h . 
we modified the initial conditions for the different dark energy sce- 
narios adapting the initial redshift for the initial conditions in the 
dark energy scenarios determined by the ratio of the linear growth 
factors D+(z), 



D + (zini) _ -D+,acdm(zacdm) 



(25) 



D+(0) D+,acdm(0) 

Additionally, the peculiar velocities of the particles are cor- 
rected according to the new redshift t o reflect a consist ent applica- 
tion of the Zel'dovich approximation dZerdovichl 1970h . 



x(t)=D+(t)H{t)VMQ) 



(26) 



Note that, unlike in previous works, here we do not use the approx- 
imation Q^lf for D+(t) as this would lead to small inaccuracies 
in some of the dark energy scenarios. Finally we also correct the 
velocities of the particles due to the changed displacement field at 
the new redshift according to 



— "ACDM 



D+(z ini )H( Zilli ) 



(27) 



Z)+,acdm(2Xcdm) #acdm(zXcdm) 
Therefore, all simulations start from the same random phases, but 
the amplitude of the initial fluctuations is rescaled to satisfy the 
constraints given by CMB. 

In Fig. [4] we show a density slice of depth equal to 1/64 of 
the box size through the whole box for each of the five models 
considered at z = 0. At first sight, we can see that the structures 
form in the same place in the different cosmologies since the initial 
phases are the same. Moreover, the differences among the models 
are small and cannot be seen with the eye; indeed, an accurate sta- 
tistical analysis is needed to understand the properties of the objects 
in the different models. 

In Fig. [5] we plot the star formation rate density (SFRD) as 
a function of redshift for all the models considered. The SFRD in 
general follows the growth of the perturbations as shown in Fig. [3] 
Using the outputs of simulations, we extract galaxy clusters 
from the cosmological boxes, using the spherical overdensity cri- 
terion to define the collapsed structures. We take as halo cen- 
tre the position of the most bound particle. Around this particle, 
we construct spherical shells of matter and stop when the total 
(i.e. dark matter plus gas plus stars) overdensity drops below 200 



ACDM 




Figure 4. Density slice of depth equal to 1/64 of the box size through the whole simulation box for the five different models at z = 



I I I I I 



ACOM 

RP 




Z 



Figure 5. SFRD as a function of redshift for the ACDM model (black), RP 
(blue), SUGRA (green), EQp (cyan), and EQn (red) cosmologies. 

times the mean (as opposed to critical) background density de- 
fined by QomPoc', the radius so defined is denoted with ifeoom 
and the mass enclosed in it as A/200m- We consider only halos 
that have M 20 om > 1-42 x 10 14 M o . We select and study ob- 
jects at three different redshifts, z = 0, z — 0.5, and z = 1. For 
the following analysis, we also calculate for each cluster selected 
in this way the radius at which the overdensity drops below 200 
(500) times the critical background density and denote it with R200 
(R&ao)- The corresponding mass is indicated as M200 (Afsoo). Just 
as a reference, the most massive object of all the simulations has 
M2oom = 3.15 x 1O 15 M0. The number of clusters at each redshift 
is different for each cosmology: for example, the sample at z — 
is made up by 563 clusters in the ACDM cosmology, 484 in RP, 
352 in SUGRA, 476 in EQp, and 431 in EQn. This fact directly 
reflects the different values of as and D+ leading to differences in 
the formation history of the halos. No morphological selection has 
been made on the sample considered, so that clusters in very differ- 
ent dynamical state are included. Nevertheless, it is useful to define 
a quantitative criterion to decide whether a cluster can be consid- 
ered relaxed or not because, in general, relaxed clusters have more 
spherical shapes, better defined centres and thus are more represen- 
tative of the self-similar behaviour of the dark m atter halos. We us e 
a simple criterion similar to that introduced in iNeto et alj (20071) : 
first of all we define x ff as the distance between the centre of the 
halo (given by the most bound particle) and the barycentre of the 
region included in i?200m ; then we define as relaxed the halos for 
which x ff < 0.07i?2oom. 



4 MASS FUNCTION 

A standard way to use galaxy clusters as cosmological probe is 
the study of their mass function. Since the total mass of these ob- 
jects is dominated by dark matter, it is a tracer of structure forma- 
tion in different cosmological models. In the top panel of Fig. [6] 
we plot the cumulative mass functions for the different cosmolo- 
gies at three different redshifts: z = 0, z — 0.5 and z = 1. This 
plot simply illustrates the number of halos per unit volume having 
a total mass greater than a given mass threshold. We can see that 
the shape and the properties of the mass functions are substantially 
the same at different redshifts (with the obvious exception of the 
maximum mass of the formed halos), with ACDM forming more 
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Table 2. Ratios between the number of clusters in the simulated volume for 
a given dark energy model with respect to Nacdm in the given M200m 
bin at different redshifts. 



M 20 0m [10 14 M Q ] z Nacdm RP SUGRA EQp EQn 



> 1.42 





563 


0.86 


0.63 


0.85 


0.77 


> 1.42 


0.5 


202 


0.81 


0.52 


0.80 


0.69 


> 1.42 


1 


45 


0.78 


0.47 


0.76 


0.64 


1.42 - 5 





507 


0.88 


0.65 


0.86 


0.78 


5-10 





45 


0.69 


0.42 


0.69 


0.67 


> 10 





11 


0.82 


0.36 


0.82 


0.64 



clusters of a given mass with compared to the other cosmologies; 
SUGRA is the cosmology which forms fewer clusters, while RP, 
EQp and EQn lie in between, with RP and EQp being the closest 
to ACDM. This fact seems to directly reflect the redshift evolution 
of the equation of state parameter w (see Fig. [2} and of the growth 
factor (see Fig. O, given the different value of as in the differ- 
ent models. Actually, for extended quintessence models, a positive 
value of the coupling £ leads to G > G, in the past, and vice versa 
for a negative f . Therefore, the linear density contrast is expected 
to be higher for EQp than for EQn. In a spherical collapse model 
like the Press-Schechter formalism JPress & Schechterj|l974T ), this 
implies a higher mass function for models with negative coupling 
(i.e. EQn) than for models with positive coupling (i.e. EQp), when 
all the other parameters are kept fixed. In our case, this effect is 
somehow mitigated by the different as used. 

In the bottom panel of Fig.|6]we plot (always at z — 0, z = 0.5 
and z = 1) the ratios between the number of clusters in a given dark 
energy model with respect to the corresponding value in ACDM. 
For each cosmology, we consider only bins in which we have more 
than one object. The same results are summarized also in Table[2] 
At * = 0, RP, EQp, EQn, and SUGRA form 86%, 85%, 77%, and 
63% the number of objects formed in ACDM, respectively. These 
numbers decrease with increasing redshift, reaching, at z = 1, 
78%, 76%, 64%, and 47% for RP, EQp, EQn, and SUGRA, re- 
spectively. This fact indicates that the differences in the formation 
history are more evident at high redshift. If we consider different 
mass bins at z = 0, we see that the differences between ACDM 
and the other models are enhanced for very massive objects, in par- 
ticular for SUGRA. 

Note that we have considered here minimally coupled models 
and scalar-tensor theories, as illustrated in Sect. [2] Couplings with 
dark matter only, where, as in equation J20b . an additional velocity- 
depend ent term is present, hav e been shown to lead to different 
results dBaldi & Pettorinol201 lh . increasing the number of massive 
clusters at high redshift. Differe nces between these sets of mod els 
have been illustrated in detail in Pettorino & Baccigalupi (2008). 

In principle, if we can count all the clusters above a given 
mass threshold, or in a given mass bin, we can try to discriminate 
between different cosmologies just using cluster number counts 
coming from cosmological surveys. From a practical point of view, 
evaluating the mass of galaxy clusters requires the assumption of 
some hypotheses on their dynamical state, and in general it is not 
an easy task to perform. So it is better to consider cluster properties 
that are directly observable (like X-ray luminosity and temperature) 
in order to distinguish among different cosmologies. We discuss 
these topics in the next two sections. 
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Figure 6. (Top panel) Cumulative mass function at 2 = 0, z = 0.5 and 
z = 1 for the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and 
EQn (red) cosmologies. For each cosmological model the mass of the most 
massive object at each redshift is marked by a cross. Error bars (shown only 
for ACDM for clarity reasons) are Poissonian errors for the cluster number 
counts. (Bottom panel) Ratios between the mass functions for different dark 
energy cosmologies and the corresponding values for ACDM at z = 0, 
z = 0.5 and z = 1. 



5 L-T RELATION 

Once we have analysed the general composition of our sample, we 
can now proceed with the study of the properties of the objects in- 
side the sample. We recall that, when considering self-similar evo- 
lution of gravitational systems, we can derive simple scaling rela- 
tions between their properties. The existence of such scaling rela- 
tions is confirmed by observations, even if in general they have a 
different shape compared to the ones predicted by self-similarity, 
indicating an important role of some non-gravitational physics in 
the evolution of these systems. We use our hydrodynamical sim- 
ulations in order to understand whether the baryon physics intro- 
duces any scale dependence that can break the self-similarity of the 
scaling relations. Since one of the aims of this work is to study 
whether there exist observable quantities that can be used to distin- 
guish among the different cosmologies considered, we start study- 
ing the X-ray L — T relation of our sample, also comparing it to 
observations to verify that the observed relation holds for our sim- 
ulated objects too. 

In order to do that, for each cluster we want to analyse we 



produce 2D maps of (5 Mpc) 2 size of the X-ray luminosity Lx 
and emission-weighted temperature T ew in the [0.5 — 2] keV soft 
band. The latter is defined by 

/ A(T)n 2 TdV 



J A(T)n 2 dV 



(28) 



where n is the gas density and A(T) is the cooling function. 

Then, for each object, we evaluate the total luminosity and the 
emission-weighted temperature in the region [0.15 — l]i?soo- We 
decide to cut the core for two reasons: first of all, despite the fact 
we use accurate physical models to describe the hydrodynamics of 
the simulations, still we do not include AGN feedback, so they are 
not optimized for the study of the central regions of the clusters; 
secondly, we have checked that cutting the core we obtain a lower 
dispersion of our data in the L — T plane. We stress that despite 
excluding the central region of the clusters in our analysis we can 
still draw robust conclusions from a cosmological point of view, 
avoiding the effects of detailed physical processes which can affect 
the inner parts. Moreover this cut is often used in observations to 
avoid problems with cool-core emission that can lead to a deviation 
from the self-similar scaling relation. Having generated luminosity 
and temperature catalogues of our sample, we can proceed with the 
analysis of the L — T relation. In the left panel of Fig.|7]we plot the 
L — T relation at different redshifts (z — 0, z — 0.5 and z — 1) 
for the ACDM cosmology. Here we correct the luminosity using 
E^ 1 = Hq/H, which is a factor c ontaining all the p redicted de- 
pendence on the cosmology (see e.g. lEttori et alj2004l) . We can see 
that there are not substantial differences at the various redshifts, 
but in general at high redshift we lack clusters in the luminosity 
region below 10 43 erg s _1 and in the temperature region below 2 
keV. This fact can be explained as a selection effect in our sam- 
ple: at high redshift, only more evolved (and thus more luminous 
and hotter) clusters are massive enough to be included in our sam- 
ple. We also provide a fit to our points, fitting the linear relation 
between the logarithms of luminosity and temperature. We find a 
slope of 1.81, which is slightly higher than the self-similar value of 
1.5 expected for the soft band considered. Finally we pl ot a collec- 
tion of observed data at different redshifts compiled bv lPratt et al.l 
(2009). The luminosities are taken exactly in the same way as we 
did, i. e. in the [0.15— 1 ] -R500 region and in the [0.5 — 2] keV band, 
while the y use spectroscop ically determined temperatures (see the 
details in IPratt et al.l 120091) . The slope of their best-fit relation is 
2.53 ± 0.16, steeper than what we found. Despite the difference in 
the slope, we can see that in the high-temperature/high-luminosity 
region where we have a sufficient number of both observed and 
simulated objects, the agreement is very good. In any case, we 
stress that a direct comparison between simulations and observa- 
tions is not the main target of this work. Here, we just want to show 
that our simulated clusters lie in a region in the L — T plane which 
is the same as the observed objects. Regarding the differences we 
find in the low-temperature/low-luminosity region, we stress that it 
is not due to overcooling in the simulations, since we are cutting 
the core; more likely, this region is populated by objects with lower 
mass, for which the detailed physical processes acting in the in- 
ner regions (e.g. AGN feedback) have im portant effects also on the 
overall properties of the clusters (see e.g. IPuchwein et al]|2008l) . 

In the right panel of Fig. [7] we plot the evolution with redshift 
of the mean luminosity and temperature in the different cosmolo- 
gies. We consider only the relaxed clusters at z = 0, z — 0.5 and 
z = 1. Then, for each cosmology, we select the ten most massive 
objects at each redshift, using M200 for this selection. Actually, at 
z = 1 for the SUGRA model we only have six relaxed clusters, 
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Figure 7. (Left panel) The X-ray L — T relation in the [0.5 — 2] keV band, evaluated in the [0.15 — 1]-Rsoo region at z = (black), z = 0.5 (cyan), and 
z = 1 (red) for the ACDM cosmology. The dashed black line is the best-fit of our simulated data. The violet squares are a collection of observed data from 
Pratt et al. (2009), while the violet line is their best-fit of the same dataset. (Right panel) Redshift evolution of the mean luminosity and temperature in the 
[0.5 — 2] keV band for the ten most massive relaxed objects in the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn (red) cosmologies. 
Circles refer to objects at z = 0, triangles to objects at z = 0.5 and squares to objects at z = 1. Dashed lines and empty symbols indicate the evolution of 
the mean luminosity and temperature evaluated inside -R500. while solid lines and filled symbols refer to the same quantities evaluated in the [0.15 — lj-Rsoo 
region. In both panels, the cosmological dependence is taken into account using the factor EJ = Ho / H which multiplies the luminosity. 



and we consider all of them. At this point, at each redshift, we eval- 
uate the mean luminosity and temperature of the selected objects 
both in the region inside -R500 and in the region [0.15 — l]i?5oo- 
We find that cutting the core results in both a lower mean luminos- 
ity and lower mean emission-weighted temperature. As a general 
trend, either including or cutting the core, both the mean luminosity 
and temperature increase with decreasing redshift, independently 
of the cosmological model. This is in somehow expected, since at 
late cosmic time the clusters are more evolved, and thus hotter and 
more luminous. The differences in the values of mean luminosity 
and temperature among the different cosmologies reflect the differ- 
ent histories experienced by objects in different dark energy envi- 
ronments, substantially following the mass function. 



6 X-RAY OBSERVABLE FUNCTIONS 

Using the same maps built to study the X-ray L — T relation, we 
can also analyse the X-ray luminosity function (XLF) and the X-ray 
temperature function (XTF) of our samples. Since the samples are 
mass selected (see Sect.[4}, only the mass functions we have shown 
before can be considered complete. XLFs and XTFs in a sense re- 
flect the mass functions, but cannot be considered complete for the 
selection effect discussed in the previous section. This means that 
at higher redshift, we are missing more and more clusters in the 
low-luminosity region of the XLF and in the low-temperature re- 
gion of the XTF. We show in Fig.[8]the cumulative XLFs and XTFs 
of our sample at z = 0. We cut the plots at 0.1 x 10 44 erg s _1 
and 1 keV in order to be as complete as possible also in the low- 
luminosity and low-temperature regions. In the left panel of Fig. [8] 
we show the cumulative luminosity function. In the middle panel of 
the same figure we plot the ratios between the number of clusters 
in a given dark energy model with respect to ACDM in every lumi- 
nosity bin. As in the case of the mass function, for each cosmology, 
we consider only bins in which we have more than one object. The 
results for three luminosity bins are also summarized in Table [3] 
In general, despite some noisy oscillations, the ratio is decreasing 



with increasing luminosity. Nevertheless, in the range between 0.5 
and 1 x 10 44 erg s _1 it increases and in three models out of four 
the number of objects is equal or even larger than in ACDM. This 
effect seems to be statistically significant in particular for RP. In 
any case, by looking only at very luminous objects, the differences 
with ACDM are significant for all models. 

In the right panel of Fig. [8] we show the same as in the left 
panel, but for the cumulative temperature function (see also Table 
[3). In this case, the decrease of the ratio with increasing tempera- 
ture is evident in all the dynamical dark energy cosmologies. Going 
from objects in the range between 1 and 3 keV to objects with tem- 
peratures higher than 3 keV, RP goes from 87% to 70%, SUGRA 
from 64% to 33%, EQp from 86% to 57%, and EQn from 78% to 
43%. 

In general, we see that the relative trend among the different 
cosmologies shown by the mass functions at z = is almost pre- 
served in the XLFs and XTFs: in a given mass, luminosity and 
temperature bin, ACDM forms more clusters than the other cos- 
mologies do (except for RP in a luminosity bin, as noted before). 
On the other hand, SUGRA is the cosmological model that forms 
fewer clusters in each bin. EQp and EQn lie in between. This find- 
ing is confirmed by the bottom panels of Fig. [8] where we show the 
ratios between the XLFs and XTFs plotted in the top panels and 
the ones recovered by applying the L — M relation at z = for 
the ACDM cosmology to the mass functions of each dark energy 
model. This is done to disentangle the differences in the XLFs and 
XTFs due to a different mass function and the ones due to baryon 
physics. The fact that the subsample considered in the right panel 
of Fig.|7]reproduces the XLF and XTF of Fig.[8]also seems to indi- 
cate that relaxed and massive objects are still a good representation 
of the whole sample. The general trend of the mass, luminosity and 
temperature functions seems to reflect the evolution with redshift 
of the dark energy equation of state parameter w = w(z), as we 
showed in Fig. [2] ACDM tends to form massive clusters earlier 
than the other cosmologies, thus giving a larger number of evolved 
(i.e. with high luminosity and temperature) objects at z = 0. The 
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Figure 8. (Left panel) The X-ray luminosity function (XLF) in the [0.5 — 2] keV band, evaluated in the [0.15 — lJ-Rsoo region at z = for the ACDM (black), 
RP (blue), SUGRA (green), EQp (cyan), and EQn (red) cosmologies. For each cosmological model the luminosity of the object with the highest luminosity is 
marked by a cross. Error bars (shown only for ACDM for clarity reasons) are Poissonian errors for the cluster number counts. In the middle panel the ratios 
between the luminosity functions for different dark energy cosmologies and the corresponding values for ACDM are shown. In the bottom panel we plot the 
ratio between the luminosity functions shown in the top panel and the ones recovered by applying the L — M relation at z = for the ACDM cosmology to 
the mass function of each dark energy model. (Right panel) The same as in the left panel, but for the X-ray temperature function (XTF). 



Table 3. Ratios between the number of clusters in the simulated volume for 
a given dark energy model with respect to N acdm m the given Lx , T ew , 
M gas500 , and Yxsoo bin - 



z = Nacdm RP SUGRA EQp EQn 
L x [lO^crgs" 1 ] 

0.1 -0.5 391 0.87 0.61 0.86 0.78 

0.5- 1 24 1.38 0.63 1.13 1.00 

> 1 13 0.62 0.23 0.46 0.46 

T ew [keV] 

1 -3 528 0.87 0.64 0.86 0.78 

> 3 30 0.70 0.33 0.57 0.43 

M gaa500 [1O 13 M ] 

1 -5 347 0.85 0.58 0.81 0.76 

> 5 14 0.50 0.21 0.43 0.50 

Yxsoo [10 13 M Q keV] 

1 -5 392 0.87 0.63 0.82 0.79 

5 - 10 55 0.98 0.49 0.96 0.96 

> 10 29 0.69 0.38 0.62 0.66 



XTF seems to better reflect the mass function, while the XLF is 
more influenced by baryonic physics, as we can clearly see from 
the behaviour of the RP cosmology. So, in principle, we can try 
to distinguish among different cosmologies by building the XTF 
of a sample of galaxy cluste rs. The problem is that if we check, 
for example, the sample from lPratt et alj d2009f) . there are very few 
clusters in the temperature range we have considered for our XTF. 
Being an X-ray selected sample, in general they have a higher tem- 
perature compared with our simulated objects, and so it is not easy 
to directly compare our results with their observational data. 

In order to check whether other proxies could better trace 



the formation history of structures, we also analysed the X-ray 
M aas5 oo and Yxsoo functions. M gas soo is defined simply by the 
mass of X-ray emitting gas contained in R500, while Yxsoo = 
M sas 5oo x Tew, where T ew is evaluated in the [0.15 — l]-Rsoo 
region. We see from Table [3] that, for M gas50Q > 5 x 1O 13 M , 
M gas 5oo is in principle a very powerful tool to distinguish between 
different cosmologies. In fact, all the models form at most 50% the 
number of objects formed by ACDM, and, since M gas soo is quite 
an easy quantity to estimate from observations, it should be pos- 
sible to rule out some models just by studying the Mg as mo func- 
tion. The quantity Yxsoo does not seem to be as good as M ga s500, 
since the differences between ACDM and the other models are less 
pronounced, and also the behaviour in the different bins is not so 
smooth. 

It is interesting to evaluate the volume that a cluster sur- 
vey must cover to be able to discriminate using the local (i.e. 
at z = 0) cluster counts among the different dark energy mod- 
els here considered. For that we assume Poissonian error bars 
and consider a 3<r level. Using the mass function with a thresh- 
old of 5 x 1O 14 M , we find that cosmological volumes larger 
than 1.6 x 10 7 (Mpc ft" 1 ) 3 are sufficient to distinguish between 
SUGRA and ACDM, while 6.4 x 10 7 (Mpc ft" 1 ) 3 are required 
for EQn and 9.1 x 10 7 (Mpc ft" 1 ) 3 are required for RP and EQp. 
Considering the XLF (with a threshold of 0.5 x 10 44 erg s" 1 ), 
larger surveys are required: in fact volumes larger than 4.3 x 10 7 , 
3.4 x 10 8 , 1 x 10 9 , and 1.3 x 10 9 (Mpc ft" 1 ) 3 are necessary to 
discriminate among SUGRA, EQn, EQp, and RP and ACDM, re- 
spectively. The situation is better when the XTF (with a threshold 
of 3 keV) is used: predictions for the ACDM model are different 
at 3(J level with respect to the ones for SUGRA, EQn, EQp, and 
RP, when volumes as large as 2.7 x 10 7 , 4.3 x 10 7 , 6.4 x 10 7 , 
and 1.7 x 10 8 (Mpc ft" 1 ) 3 are considered, respectively. This fact 
confirms the importance of XLF/XTF in tracing the number counts 
in a given cosmology, and again that the XTF is a better quan- 
tity to be used in that kind of studies, if compared to the XLF. 
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We recall that we are not considering any selection function on 
XLF/XTF, being a proper treatment of any observational approach 
in defining an XLF/XTF beyond the purpose of the present work. 
If we move to z = 1, using the mass function with a threshold 
of 1.42 x 1O 14 M , 2.7 x 10 7 (Mpc ft" 1 ) 3 are still sufficient to 
distinguish between SUGRA and ACDM, while EQn, EQp, and 
RP need 6.4 x 10 7 (Mpc /i" 1 ) 3 , 1.7 x 10 8 (Mpc ft" 1 ) 3 , and 
2.2 x 10 8 (Mpc h' 1 ) 3 to be distinguished from ACDM, respec- 
tively. Larger cosmological boxes and larger observational samples 
with higher resolution and sensitivity (i. e. lower flux limit), such 
as, e.g., the one expected with eROSITA JPredehl et alj|2007h and 
WFXT l lGiacconi et aLl l2009l) . can provide better answers to the 
question. 
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7 THE BARYON FRACTION 

In this section we focus on the baryon fraction ft ar — f s tar + fgas 
of our simulated galaxy clusters, where f 3 t ar = M s tar/Mtot and 



fga 



Mgag/Mtot- Since we are considering galaxy clusters in 



a cosmological context, it is better to re-express the star fraction 
fstar, the gas fraction / JM , and the total baryon fraction ft ar in 
units of the cosmic baryon fraction fiob/^Om = 0.164 adopted in 
these simulations. To do this we introduce the quantities 



fstar 



fgas 



fbar 



Qot I fit) II 



n 



■n 



,(29) 



and indicate them as star, gas and baryon depletion parameter, re- 
spectively. In this section we analyse the dependence of these quan- 
tities on mass, redshift and distance from the centre of the object 
considered, as well as on the underlying cosmology. 

In Fig.[9]we plot the distribution of b oar , b gas and b s tar eval- 
uated at 7?2oo for the whole sample at z — 0, z — 0.5, and z = 1 
in order to check the spread of the values for the single objects. 
We see that at z = there is a substantial overlapping among the 
different cosmologies, indicating that evolved objects have almost 
the same distribution whatever the underlying cosmological back- 
ground is. The same is true looking at z = 0.5 and 2 = 1. We can 
note a decrease of b ga s going from z = 1 to z = 0, not completely 
compensated by an increase of b s tar- The net effect is a decrease 
of bbar going from z — 1 to z = 0. 

In Table [4] we summarize the mean value of b s tar, b ga s, and 
bbar evaluated at -R2500, R500, and -R200 for all the objects in the 
different cosmological models considered, at z — 0, z — 0.5, and 
a = l. We see that, on the one hand, for any cosmological model, 
at any redshift, b s tar is a decreasing function of radius, going from 
/fesoo to 7?200- On the other hand, b ga3 is an increasing function of 
radius. As a whole, b aar is slightly decreasing with radius. Fixing 
the radius, either -R2500, R&oo, or R200, b s tar increases going from 
z = 1 to z = 0, while b ga s decreases. As we already noted from 
Fig. [9] btar is slightly decreasing going from z = 1 to z = 0. 
These trends are general, and they hold for all the cosmological 
models considered. 

In the left panel of Fig.[l0]we plot, for each cosmology, the 
ratio between the mean values of bb ar , b gas and b s tar evaluated at 
R200 in four different mass ranges at z = and the mean value 
of bZr m , b™ r s m and 6™ t °™ for clusters having A/200 > 1O 15 M 
in the ACDM cosmology at z = (i.e. 0.921, 0.757 and 0.165 
respectively). We have considered four mass ranges: A/200 < 
10 14 M Q , 1O 14 M < A/200 < 5 x 1O 14 M , 5 x 10 14 M Q < 
A/200 < 1O 15 M , and M 20 o > 1O 15 M . We have evaluated the 
quantities at R200 instead of iisoo as in Sect.|5]because this radius 
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Figure 9. Distribution of fcf,ar> b ga s, and b s tar evaluated at -R200 f° r the 
ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn (red) 
cosmologies at z = (top panel), z = 0.5 (middle panel), and z = 1 
(bottom panel). 



is representative of the cluster as a whole, including the external 
regions, and indeed we want to check whether, in different cos- 
mologies, these objects are a fair representation of the underlying 
background. The first thing we notice is that, in every mass bin, the 
values of bbar, b g as and b s tar are similar, within error bars, among 
the different cosmologies. So we can refer to a single cosmology 
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Figure 10. The evolution of stellar, gas and baryon depletion parameter evaluated at -R200 with mass at 2 = (left panel) and 2 = 1 (right panel) for 
the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), and EQn (red) cosmologies. Crosses, triangles, and squares indicate b s tar, b gas , and bf, ar 
respectively. The depletion parameters are expressed in units of b^°£ m , bg° B m and b™°™, the mean values for clusters with M200 ^ 10 15 Mq in the 
ACDM cosmology at 2 = 0. Error bars are r.m.s. of the mean distribution. 



Table 4. Mean values of b s tar, b gaB , and b^ ar evaluated at -R2500. ^500. an d ^200 f° r all the objects in the different cosmological models considered, at 
2 = 0, 2 = 0.5, and 2 = 1. N indicates the number of objects in the given model at the given redshift. Numbers in brackets are la errors on the mean. 



Model 


z 


N 


bstar2500 


bgas2500 


bbar2S00 


bstar500 


&gas500 


&bar500 


bstar200 


bgas200 


bbar200 


ACDM 





563 


0.508 


0.535 


1.043 


0.269 


0.680 


0.948 


0.207 


0.731 


0.937 








(0.072) 


(0.067) 


(0.065) 


(0.032) 


(0.037) 


(0.036) 


(0.021) 
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(e.g. ACDM) in order to study the mass dependence of these quan- 
tities. We see that bbar is almost constant, independently of mass. 
On the one hand, b ga s shows a slight positive trend, of the order 
of 5%, going from low-mass to high-mass systems, but still com- 
patible with a constant value within the error bars. On the other 
hand, b ata r shows a decrease up to 30% going from low-mass to 
high-mass clusters, not compatible with a constant value. In the 
right panel of Fig. [10] we plot, for each cosmology, the ratio be- 
tween bb ar , b gas and b 3ta r evaluated at -R200 at z = 1 and the 
mean values for ACDM at 2 = already used in the left panel. 
In this case we do not consider different mass ranges, since at this 
redshift the cluster abundance starts to be low and all the objects 
have 1O 14 M < M 20 o < 5 x 10 14 M Q . A gain, the cosmologies 
are completely equivalent within error bars. 

Here we stress again that our simulati ons do not follow AGN 
feedback. It is known from literature (e.g. iPuchwein et alj|2008l) 
that the effect of this feedback is mass dependent, leading to a low- 
ering in the baryon fraction in groups and low-mass clusters, with- 
out affecting significantly high-mass clusters. 

We find in general a co nstant baryon fractio n with respect to 
the mass. Some authors (e.e. iGiodini et alj|2009h claim that in ob- 
served objects the total baryon fraction shows an increase with in- 
creasing mass. This difference with respect to our results could be 
due to the fact that some relevant physical processes are still not 
included in our cosmological simulations. Such processes may be 
able to affect the global properties of low-mass systems without 
changing the high-mass clusters. Not including them in the simu- 
lations does not permit to us to fully compare our results with ob- 
servations. In particular, we note an overabundance of stars (which 
obviously influences the total baryon fraction) in low-mass objects. 

Combining the right and left panels of Fig. [TO] we can study 
the evolution with redshift of bbar, b gas and b sta r. Since the differ- 
ences among various cosmologies at the same redshift are quite 
small, we rely on our reference ACDM model for the analysis 
of redshift evolution. For clusters with 10 14 Mq M200 < 
5 x 1O 14 M0, the mean value of bb ar is almost constant, with a 
slight decrease of about 2%, going from z = 1 to z = 0. In par- 
ticular, bg as decreases of less than 10%, while the increase of bstar 
is of the order of 25%. A decrease of the baryon fraction with de- 
creasing redshift was already found in other simulations (see e.g. 
lEttori et aljbood) . and a possible explanation is that at high red- 
shift the radius at which the baryons accrete is smaller than at low 
redshift, and so a greater number of baryons can fall in the cluster 
potential well. 

Finally, we study the star, gas, and baryon depletion param- 
eters as a function of the distance from the centre of the cluster, 
defined as the position of the most bound particle. For each cos- 
mology we select, as in Sect. [5] the ten most massive (in M200) 
relaxed halos and generate the radial profile of the object obtained 
by stacking them. We do this at z = 0, z — 0.5 and z = 1. We 
recall that, at z = 1, for SUGRA we only have six objects. The 
resulting profiles, expressed in units of the cosmic baryon frac- 
tion flob/^Om = 0.164, are shown in Fig. QT] At z = 0, in the 
outer regions near -R200, the five cosmologies are completely equiv- 
alent, with bbar having almost the cosmological value, while look- 
ing toward the centre some differences can be seen. This fact means 
that, as a whole, evolved relaxed clusters contain the same amount 
of baryonic matter, independently of the underlying cosmological 
model, but that the matter can be redistributed inside them accord- 
ing to their formation history. This fact is confirmed by looking at 
z = 0.5 and in particular at z = 1, where the differences among the 
models are clear even in the outer regions, indicating a sort of self- 
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Figure 11. The evolution of stellar, gas and baryon depletion parameter 
with radius at z = (top panel), z = 0.5 (middle panel), and 2 = 1 (bot- 
tom panel) for an object obtained by stacking the ten most massive relaxed 
clusters in the ACDM (black), RP (blue), SUGRA (green), EQp (cyan), 
and EQn (red) cosmologies. Crosses, triangles, and squares indicate b B tar, 
b g as, and bb ar respectively. The depletion parameters are expressed in units 
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the mean distribution. 
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regulating mechanism that leads to the same objects at z = even 
if they can be very different at higher redshifts. Again, the same 
features appear both in the mean values of the whole sample and 
in more relaxed and massive objects, indicating that the latter are a 
fair representation of the clusters in a given cosmological model. 

As a general rule for the radial profiles, it is confirmed the 
well known relative trend of the radial profile of gas and stars 
components, being the former increasing with radius and the lat- 
ter decreasing, giving a total baryon fraction almost constant (but 
slightly decreasing) beyond 0.5i?2oo- Then we note that the total 
baryon fraction at z = 1 is higher compared to z = 0, in partic- 
ular in the inner regions of clusters. The effect is mainly due to a 
higher star fraction in the inner regions of clusters at z = 1. An- 
other quite evident feature is that the radius at which the gas starts 
to dominate over the stars increases with increasing redshift. The 
explanation is that, as we have just seen, the gas fraction profile is 
almost independent of redshift, while the star fraction at a given ra- 
dius increases with redshift, and so at higher redshift it remains the 
dominant baryonic component also at larger radii. 



8 CONCLUSIONS 

In this paper we have analysed the general properties of a sam- 
ple of galaxy clusters extracted from hydrodynamical simulations 
of different cosmological models with dynamical dark energy. We 
simulate a cosmological box of size (300 Mpc h' 1 ) 3 , resolved 
with (768) 3 dark matter particles and the same amount of gas par- 
ticles. The reference cosmology is a concordance ACDM model 
normalized to WMAP3 data. The others are two quintessence mod- 
els, one with a RP and the other with a SUGRA potential, and two 
extended quintessence models, with a positive and a negative cou- 
pling between quintessence and gravity, indicated as EQp and EQn, 
respectively. Since all models are normalized to CMB data, they 
have different erg, and thus different structure formation histories. 
We focus on various properties of the considered objects, in partic- 
ular the mass function, the X-ray L — T relation, the X-ray lumi- 
nosity and temperature functions (XLF and XTF respectively) and 
finally the baryonic fraction in terms of the depletion parameters 
bstar, bgas and bt, ar defined in equation \29\ . We select and study 
objects at three different redshifts, z — 0, z — 0.5, and z = 1, 
with M200W1 ^ 1-42 x 1O 14 M0. We also define a criterion to dis- 
tinguish between relaxed and unrelaxed clusters. From our analysis 
we draw the following conclusions. 

• Mass function: at 2 = the total mass function evaluated at 
R200m shows different behaviours in the different cosmologies, in 
particular in the normalization. The ACDM model tends to form 
more clusters of a given mass with compared to the other cosmolo- 
gies; SUGRA is the cosmology with the smallest abundance, while 
RP, EQp and EQn lie in between, with RP and EQp being the clos- 
est to ACDM. This fact directly reflects the redshift evolution of 
the equation of state parameter w and of the growth factor, given 
the different assumed value of as in the various models. Actually, 
for extended quintessence models, a positive value of the coupling 
leads to an higher linear density contrast, and vice versa for a neg- 
ative coupling. This would imply a higher mass function for mod- 
els with negative coupling (i.e. EQn) than for models with positive 
coupling (i.e. EQp), keeping fixed all remaining parameters. In our 
case, this effect is somehow mitigated by the different erg used. 
This trend is preserved also at z — 0.5 and 2 = 1. The differ- 
ences among the models are more pronounced in the high-mass tail 
of the distribution. This is expected, because very massive objects 



form later and are representative of the different structure forma- 
tion time scale of the considered cosmologies. We note here that 
our results are different fr om what has been foun d in the case of 
coupling with dark matter JBaldi & Pettorinol201 lh . where there is 
an enhancement in the number counts of massive objects. 

• L — T relation: we compare the L — T relation of our sim- 
ulated objects in t he ACDM refer ence models with a collection of 
observed objects jPratt et aT]|2009h . Despite the differences in the 
slope of the relation in the two cases (1.81 for our simulated objects 
vs 2.53 ± 0. 16 for their observed ones), we find that there is a good 
agreement in the high-temperature high-luminosity region where 
X-ray selected observed objects are found. The discrepancy in the 
low-temperature low-luminosity region is not worrying, because 
low-mass systems are globally more affected by phy sical mecha- 
nisms not yet fully understood and reproduced (e.g. [Borgani et al.l 
2004), acting in the core. We also study the evolution with redshift 
of the L — T relation for the ten most massive relaxed objects in 
each cosmology, both keeping and cutting the core. We find that 
cutting the core results in both a lower mean luminosity and lower 
mean emission-weighted temperature. In general, both the mean 
luminosity and temperature increase with decreasing redshift, in- 
dependently of the cosmological model, because they trace the hi- 
erarchical growth of structures. 

• X-ray observable functions: the relative behaviour observed in 
the mass functions is also qualitatively reproduced by the XLFs and 
XTFs evaluated in the [0.5 - 2] keV band in the [0.15 - l]R 50 o 
region, with few exceptions. In particular, in the range of luminos- 
ity around 0.5 x 10 43 erg s _1 RP tends to form 10% more clusters 
than ACDM. We also check the X-ray M 90S soo and Yxsoo func- 
tions as proxies for the mass function. We conclude that all the X- 
ray observable functions are more or less equivalent, with T ew and 
M gas 5oo being slightly more stable than Lx and Yxsoo, in trac- 
ing the mass function and thus disentangle the growth of structures 
among different dark energy models. For each dark energy model 
we evaluate the volumes that a cluster survey must cover in order to 
be able to distinguish it from the concordance ACDM model, using 
the mass function, the XLF, and the XTF. 

• Baryon fraction: the analysis of the b s t ar , b gaa , and bbar de- 
pendence on mass, redshift and distance from the cluster centre 
shows that there is no significant difference among the five cos- 
mologies considered, if we limit ourselves to the values at R200 and 
at z = 0. Therefore, at these conditions, bt ar (and so the baryon 
fraction fb ar ) can be safely used as a cosmological proxy to derive 
the value of other cosmological parameters. In addition, we do not 
find any clear positive trend of the total baryon fraction with mass, 
while we see a positive trend (of the order of 5%) of the gas fraction 
and a negative trend (of the order of 30%) of the star fraction going 
from low-mass to high-mass systems. Considering observations of 
real objects, in spite of finding the s ame trend for the g as and star 
fraction as we do, other authors (e.g. lGiodini et al . 2009) claim that 
the total baryon fraction is increasing with increasing mass. Ac- 
tually, for all the cosmological models here considered, we find a 
slight decrease in the total baryon fraction with increasing mass. 
Still, we have to recall that, despite the hydrodynamical treatments 
in the simulations is based on sophisticated physical models, we do 
not include AGN feedback in our simu lations. It is known from lit- 
erature (see e.g. IPuchwein et al] |2008) that AGN feedback is mass 
dependent, in the sense that it globally affects more low-mass sys- 
tems than high-mass systems. The net effect is the lowering of the 
total baryon fraction in low-mass objects while not affecting more 
massive clusters. Finally, we find a slight decrease (at most 5%) of 
the baryon fraction going from high to low redshift. A similar trend 
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was already noted bv lEttoriet all [2006) and a possible explanation 
is that at high redshift the radius at which the baryons accrete is 
smaller than at low redshift, and so a greater number of baryons 
can fall in the cluster potential well. 

In the end, we can conclude that in models with dynamical 
dark energy, the evolving cosmological background leads to differ- 
ent star formation rates and different formation histories of galaxy 
clusters, but the baryon physics is not affected in a relevant way. In- 
deed, evolved and relaxed clusters, if studied in regions sufficiently 
far from the centre, reveal to be very similar despite the different 
dark energy models considered. So, in conclusion, galaxy clusters 
can effectively be used as a probe to distinguish among different 
dark energy models. 
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